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ABSTRACT 

The effect of suction on controlling the three-dimensional secondary insta- 
bility is investigated for a boundary layer with pressure gradient in the pres- 
ence of small but finite amplitude Tollmien-Schlichting wave. The focus is on 
principal parametric resonance responsible for strong growth of subharmonics 
in low disturbance environment. Calculations are presented for the effect of 
suction on the onset and amplification of the secondary instability in Blasius 
and Falkner-Skan flows, as well as its effect on controlling the production of 
the vortical structure. 
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1. INTRODUCTION 

Laminarization of the flow by suction, and subsequent viscous drag reduction, is the prin- 
cipal and most effective means used for laminar flow control (LFC) [1]. Excessive suction 
increases suction drag and reduces the overall efficiency of an aircraft with LFC. Moreover, 
the over-thinned boundary layer is over-sensitive to surface irregularities. Hence, it is neces- 
sary to keep the boundary layer laminar with the least possible suction. For this purpose, the 
stability characteristics of the flow need to be accurately calculated, and the laminar-turbulent 
transition process must be well understood. Moreover, the effect of suction control on some of 
the early stages to transition need to be assessed. 

The nonlinear process of the flow once initiated becomes violent and lead rapidly to tran- 
sition, and it is then extremely difficult to control the flow field. Hence, suction control in 
LFC systems is preferred to be in the linear range. For these systems to be efficient, some dis- 
turbance growth should be allowed, and linear primary stability theory in conjunction with the 
e" method were relied on for such predictions [2]. A linear growth of a primary two- 
dimensional (2D) Tollmien-Schlichting (TS) wave may parametrically excite a linear secondary 
growth with three-dimensional (3D) character. This secondary instability may not lead to tran- 
sition by itself, but as it grows, it interacts with both the mean flow and the TS wave leading 
rapidly to transition. The effect of suction on the primary TS wave is well established and 
known to be drastic [3]. While, in this paper, we are concerned with the effect of suction on 
controlling the onset and growth of the secondary 3D instability. 
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Following the growth of the primary 2D wave, the flow takes on an increasingly three- 
dimensional (3D) behavior [4] with explosive growth rates. In this early stage of transition, 
the 3D phenomenon is characterized by a strong secondary instability with respect to 3D dis- 
turbances in the presence of small but finite amplitude primary 2D waves [5,6]. Primary 3D 
disturbances might be stable or very slowly growing in the absence of the 2D waves. This 
secondary instability has been recognized experimentally (Refs. 7-12) and observed numerically 
(Refs. 13-18) in the uncontrolled boundary layer (Refs. 7-15) as well as the controlled boun- 
dary layer (Refs. 16-18). 

Two major types of breakdown have been identified, a fundamental breakdown (K-type) 
and a subharmonic breakdown (H/C-type). The experiments indicate that the subharmonic 
breakdown occurs when the amplitude of the primary TS wave is low or moderate, while the 
fundamental breakdown occurs for higher amplitudes. However, one type or a mixture of both 
will appear depending largely on the spectrum of the background disturbances [19]. The linear 
secondary instability theory formulated by Herbert [20] can predict the increasingly 3D 
behavior with large growth rates that occur in both the fundamental and subharmonic types of 
breakdown, while Craik’s resonant triad model [21] predicts an instability of the subharmonic 
type. Craik’s mechanism (referred to as C-type) is thought to dominate at low amplitude of 
the TS wave, while Herbert’s subharmonic mechanism (referred to as H-type) reflects the situa- 
tion at moderate amplitudes of the TS wave. A recent review on the subject of secondary ins- 
tability has been provided by Herbert [22] and Bayly et al. [23]. 

The 3D subharmonic instability which characterizes the road to transition initiated from a 
low disturbance background appears to be more realistic in flight applications. Hence, in this 
paper we investigate the development of the subharmonic secondary instability in a boundary 
layer with pressure gradients controlled by suction. Our objective is to evaluate the effect of 
suction control on this early stage leading to transition. Several questions need to be answered, 
Does suction delay the onset of the secondary instability ? How sensitive is the growth of the 
secondary instability to the intensity of suction ? What is the effect of the initial amplitude of 
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the primary wave on this sensitivity ? Does the effectiveness of suction as a method for delay- 
ing transition depend on where it is applied, or on its intensity ? Finally, if the boundary layer 
would be kept laminar with the least possible suction, then should one allow for a limited 
growth of the secondary disturbance, or should one increase suction to fully stabilize the secon- 
dary disturbance ? The focus of this paper is to answer these questions as well as to shed 
some light on the mechanism by which the suction control the production of the vortical struc- 
ture of the secondary instability. 

2. ANALYSIS 

2.1. The Mean Flow 

We consider a 2D boundary-layer flow of an incompressible fluid with inviscid flow field 
given by U=U(x) and distributed suction given by v=v w (x) at the wall, where x is the stream- 
wise direction and y is the normal direction. The flow is governed by the nonsimilar 
boundary-layer equation 

fim+ffim+kl-fv HY m C 1 ) 

with boundary conditions 

Aft. 0)=0, /(t0>f2^ 4 (^0)=0, (2) 

given in the Gortler variables 

d%=U t dx, dr\=(U.lJ2£,)dy (3) 

Then the velocity components u and v in terms of the new variables are 

«=t/,/ n (4) 

v =v w -(U, /V%) [/ +24/y7i(M)/T,] (5) 

where the suction and pressure gradient functions r and (3 respectively are defined as 




( 6 ) 
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&(£M WU e 2 )(.dU e ldx) (7) 

Note that negative r indicates suction, while negative 0 indicates unfavorable pressure gra- 
dient. If both the suction and pressure gradient functions are constant (equal to r 0 and p 0 
respectively), then f(&,T\) is a function of ti only and we have a similar boundary-layer profile 
governed by 

frm+ffrm+Ml-fZ Wnn = 0 (8) 

/(0)=/ n (0) = o, y r, — >i (9) 

where the condition r<j=constant demands that v w be proportional to U e N 2£. On a flat plate, v w is 
proportional to 1 /Vx . 

2.2. The Primary Instability 

We consider the primary instability of the calculated mean flow with respect to 2D quasi- 
parallel spatially growing TS waves. Dimensionless quantities are introduced to the governing 
incompressible Navier-Stokes equations by using the reference velocity U, and the reference 
length L=(y,xlU t ) m so that Reynolds number is given by R-(U e x/v e ) m where * measures the 
distance from the leading edge, and v is the fluid kinematic viscosity. The primary TS wave is 
assumed to take the traveling form 

q(x,y,t)=A [<7i00 exp(z ja r dx-i (Of }+cc. ] + 0{A 2 ) (10) 

where 

A = A(x) = A 0 exp(-Ja, dx), (11) 

and qi stands for the velocities u, and v l5 and the pressure p x of the primary wave, and cc 
denotes complex conjugate. The spatial stability analysis is chosen for being more appropriate 
for this study. Hence, the wavenumber is complex and given by a^a r +z a, and the disturbance 
frequency co is real. The linearized incompressible Navier-Stokes equations reduce to a 
fourth-order system of ordinary differential equations with homogeneous boundary conditions 
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in U\ and v x . These equations may be combined to yield the well known Orr-Sommerfeld 
equation. The eigensolutions of the primary wave are normalized such that the amplitude A 
measures directly the maximum root mean square value of the disturbance velocity in the flow 
direction. For the purpose of comparison we follow Herbert [20] and let 

max2[«i()')] 2 = 1 (12) 

The linear stability theory of the primary wave provides a for a given co and R . Then the 
integration of the growth rate -a, gives the amplification factor In (AM 0 ), where A 0 is an arbi- 
trary initial amplitude at the onset of the primary instability. 

2.3. The Secondary Instability 

We consider a basic state q(x,y,t) which consists of a 2D quasi-parallel boundary layer 
with pressure gradient and suction modulated by a periodic component of the linear primary 
instability problem, 

q(,x,y,t) = q<ixy)+A q^xy,t) (13) 

To study the linear 3D instability of the basic state q(x,y,t), we superpose a small 
unsteady disturbance on each flow quantity of the basic state, that is 

q(x,yj,t) = q(x,y,t) + B q 2 (x,y,z,t ) (14) 

where q 2 is a secondary disturbance eigenfunction that stands for velocities u 2 , v 2 , w 2 , and pres- 
sure p 2 , they are normalized such that the amplitude B measures the maximum root mean 
square value of u 2 . The amplitude B is assumed small compared to the primary amplitude A, 
such that the primary instability will influence the secondary instability but not vice versa. 

Herbert [22] has pointed out that the 3D secondary instability occurs at small amplitudes 
of the primary wave where the nonlinear distortion is weak. This instability is of vortical 
nature and originates from a strong mechanism of combined tilting and stretching of the vor- 
tices [5], leading to large growth rates when compared with those for the primary wave. In 
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view of that, we neglect in this analysis the nonlinear distortion of the eigensolutions at 
finite amplitude of the primary wave. Also, the growth of the primary TS wave which occurs 
on viscous scales can be considered weak, and the variation of the primary amplitude A can be 
assumed locally constant. 

Equation (14) is substituted into the dimensionless Navier-Stokes equations. The basic 
state is subtracted, and the resulting equations are linearized with respect to the secondary 
amplitude B. We end up with four coupled partial differential equations for the secondary 3D 
instability. The coefficients of these equations are function of the basic state, they are indepen- 
dent of the coordinate z, and periodic in * and t. Hence, the z -variation can be separated, and 
Floquet theory of differential equations with periodic coefficients can be applied to give a solu- 
tion in the form 

q2(.x,y,z,t) = e' IXH ”e i ^ <i>(x,y,0 (15) 

where p is a real spanwise wavenumber of the secondary disturbance. Y=Yr+»Y; and <j=o r +i o, 
are two complex characteristic exponents, and <j >(x,y,t) is a periodic function of (x-m/a), the 
same as the period of the basic state. We express $ in terms of Fourier series to obtain the fol- 
lowing expression for q 2 

q&V *4) ****«*€'+ £<t >„(y) (16) 

Equation (16) represent the general Floquet form for the eigenmodes of a periodic basic 
state, where both the fundamental and subharmonic modes are special cases. Only two of the 
four real exponents y r ,y, , o r , and a,, are determined as a solution of the eigenvalue problem, 
others must be given. For the purpose of our study of the spatial instability of subharmonic 
modes, we let y r represent the growth rate of the secondary instability, a r =0 (no temporal 
growth), a,=-co/2 for pure subharmonic mode, and let y. represent the shift in the streamwise 
wavenumber of the secondary disturbance with respect to the primary one. y ,=0 means that the 
secondary disturbance is perfectly synchronized with the basic state. 
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We use the lowest truncation of the Fourier series for the subharmonic mode («<1) in Eq 
(16), to derive four coupled ordinary differential equations for u 2 , v 2 , w 2 , and p 2 . These equa- 
tions when supplemented with homogeneous boundary conditions, they constitute an eigenvalue 
problem 

y=r(a,MM) (17) 

for a given boundary-layer velocity profile with pressure gradient and suction. The eigenvalue 
problem (17) of the secondary instability provides y for a given p and R. Then the integration 
of the growth rate y r gives the amplification factor In (B/B 0 ), where B 0 is an arbitrary initial 
amplitude at the onset of the secondary instability. 

3. NUMERICAL TREATMENT 

Similar boundary layer profiles were calculated using Eqs (8)-(9). These equations were 
numerically integrated by using a shooting technique with fourth-order Runge-Kutta and 
Adams-Moulton integrator. In cases of continuous suction, where v w = constant, similar solu- 
tions do not exist, and Eqs (l)-(2) were numerically integrated by using a step by step pro- 
cedure in the streamwise direction. A three-point implicit finite difference technique was used 
to reduce them to a set of simultaneous tridiagonal equations. These equations were linearized 
and then solved using the algorithm of Thomas. The method of solution closely parallel that 
of Price and Harris [24]. 

The primary instability which modulate the 2D boundary layer is governed by a fourth- 
order system of equations. While the eigenvalue problem (17) describing the secondary distur- 
bance, is governed by a sixth-order system of equations. Both can be written as complex sys- 
tem of linear first-order ordinary differential equations in the form 

Primary 

DZin-YjanmZi m = 0, n=l,2..4 

m=l 


(18) 


where Z u are defined by 


Secondary 

Zll - Wi, Zi 2 -Du i t Z 13 -v 1% Z\4-p\ 
6 6 



OZm - J^ b nm^ 2 m =A ^c nm Z 2m , n = 1,2..6 

m= 1 m=l 

(21) 


Z<i\ — Z 22 = Z 25 = 0 at y — 0 

(22) 

where are defined by 

^ 2 i> ^23» ^25 — >0 as y —>° © 

^21 ~ U 2> Z 22 — Du 2 , Z 2 3 — V 2, 
2 24 = P 2 » ^25 = ^2, Z 25 = Dw 2 

(23) 


and the overbar denotes a complex conjugate. The a m are the elements of 4x4 variable 
coefficient matrix of the primary eigenvalue problem, while and c m are the elements of 
6x6 variable coefficient matrices of the secondary eigenvalue problem. These elements are 
given in the Appendix. Note that in case of no modulation of the mean flow by a primary 
wave ( i.e . A=0), then the system of equations (21)-(23) will govern the stability of a primary 
subharmonic 3D wave. 

Both the primary and secondary system of equations are numerically integrated as initial 
value problem using a freestream solution as initial condition. For the secondary eigenvalue 
problem, we assume that the amplitude of the primary vanishes far in the freestream at y > y e 
(e denotes the edge of the boundary layer). Then the system (21) will have constant 
coefficients and can be solved analytically producing three linearly independent exponentially 
decaying solutions to conform with the boundary condition (23). With the freestream solution 
as initial condition, Eqs (21) are integrated from y=y, to y=0 at the wall, using a variable step- 
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size algorithm [25], based on the Runge-Kutta-Fehlburg fifth-order formulas. The solution is 
orthonormalized at a preselected set of points using a modified Gram-Schmidt procedure. A 
Newton-Raphson technique is used to iterate on the eigenvalue to satisfy the last wall boundary 
condition within a specified accuracy of 0(l(T 5 ). 

4. RESULTS AND DISCUSSION 

For the Blasius and Falkner-Skan profiles, our results were checked with that of Bertolotti 
[26], and found to be in full agreement. 

All results reported here are for the nondimensional frequency F=10 6 cd //?=60, that remains 
fixed as a wave with fixed physical frequency travels downstream. For Blasius flow, a primary 
TS wave with this frequency grows between R,= 554 and R„~ 1052 (first and second neutral sta- 
bility points), reaching a maximum amplification factor of A/A 0 =42. Also in this region, primary 
3D subharmonic waves are subject to amplification for a broad band of spanwise 
wavenumbers, but the time and length scales of these instabilities are too small to compare 
with experimentally observed transition. A strong growth of secondary 3D subharmonics can 
be due to parametric excitation by the finite amplitude primary TS wave [20]. 

4.1. Flat Plate Boundary Layer Controlled by Suction 

At rt=l050 ( almost at R,i for Blasius flow), Fig (1) shows the growth rates of the secon- 
dary 3D subharmonic disturbances as function of the spanwise wavenumber parameter b and 
for various amplitudes A of the primary wave. The parameter b defined as fc=10 3 p/ R 
represents a fixed physical spanwise wavenumber for a wave traveling downstream. Figure (1) 
compares the results of a flat plate boundary layer with suction (rv=-.l) and with no suction. 
At fixed r 0 . the figure indicates a stabilizing effect on the secondary instability as the primary 
amplitude A decreases. When the amplitude is very small, considerable growth rates exist in a 
small band of spanwise wavenumbers. The maximum growth occurs for a wave with slightly 
lower spanwise wavenumber as the amplitude A decreases. This shift to lower spanwise 
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wavenumber is more noticeable with the increase of suction parameter. As Reynolds number 
increases, the growth rates of the secondary disturbances increase at fixed F and A . 

The application of active control by suction in LFC systems restrains both the mean-flow 
profile and the primary wave. We find that the prime influence of suction on the secondary 
instability is due to the appreciable decay of the amplitude of the primary wave. This conclu- 
sion is explained next with the use of Fig. (2). In Blasius flow, a primary wave with initial 
amplitude A 0 =0.00025 at R, grows such that at R =1050 downstream its amplitude reaches a 
value of 0.01. Figure (2), curve a, shows a wide band of 3D subharmonics amplifying at that 
location with the most unstable at b = 0.17. The influence of a modified mean flow can be 
demonstrated by applying suction at R = 1050 with fixed A =0.01. This decreases the growth rate 
of the subharmonic and slightly limit the band of amplified spanwise wavenumbers (curves b 
and c), while the most unstable wavenumber is not affected by suction. When suction is 
applied ( r 0 = -0.05) starting upstream at R,, the mean flow is modified and the growth of the 
primary wave is slowed down, such that at R = 1050 downstream the amplitude of the primary 

i 

wave reaches a value of 0.0015. Calculations, curve d, indicate a strong stabilizing effect on 
the secondary instability with the most unstable 3D subharmonic occurring at a lower spanwise 
wavenumber. 

Previous results may be viewed as local, they only reflect the stability characteristics of 
the secondary subharmonic disturbance at a fixed Reynolds number. To model the experiment 
and evaluate the overall effect of suction on the onset, growth rate, and amplification factor of 
the secondary 3D subharmonic, we should combine the effect of increasing the amplitude A of 
the primary TS wave as well as increasing R as the disturbance moves downstream. In Fig 
(3), we do that and show the variation with R of the growth rates of the secondary subhar- 
monic wave. In these calculations, the initial amplitude of the primary wave is assumed 
Ao=.001, and the spanwise wavenumber b=. 15 which is an average value of the most unstable 
wavenumbers for the parameters under investigation. The growth rates of the corresponding 
primary waves are also included in the figure for comparison. 
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The micro-events which lead to the secondary 3D instability go as follows. Initially, the 
primary instability sets in at R t on the primary neutral curve, and -a, starts to grow. At a cer- 
tain value of the primary amplitude, a secondary 3D mode is induced and sets in at Rj on a 
secondary neutral curve. y r starts to grow strongly due to the increase in both A and R . Ulti- 
mately, -a; begins to decay, and the primary amplitude reaches a maximum when -a,=0. 
Shortly downstream, y r reaches a maximum and then starts to decay with the secondary ampli- 
tude reaching its maximum when y r =0. The overall effect of small suction rates is to delay the 
onset of the secondary instability (occurs later downstream), and to decrease significantly its 
growth rate. Calculated amplification factors for the secondary disturbance (not shown) indi- 
cate 71 % reduction due to the increase of r 0 from 0 to -.05. Increasing r 0 to -.1 dampened 
completely the subharmonic secondary disturbance in spite of considerable growth shown by 
the primary wave. We know from the experiment [9] that the secondary instability depends on 
the primary wave amplitude as well as on the wave fetch. In a fixed disturbance environment, 
small suction rates may not affect the initial amplitude of the primary wave but certainly will 
influence its growth. Our calculations in Fig (3) indicate a delay in the onset of the secondary 
instability due to suction. This delay is accompanied by a slight decrease in the primary thres- 
hold amplitude (equal to .0029 at r^O and .0024 at r^-.OS). However, at r<p-.i, the primary 
amplitude reaches a maximum of only .0014 which is apparently below the value needed to 
induce a secondary subharmonic instability. Notice that the initial amplitude of the primary 
wave is fixed in all cases and equal to .001. 

Figure (4), a case of suction rate ro=-.05, shows a primary instability that sets in at R =650. 
The onset of the subharmonic instability occurs at R = 850 when Ao=. 001, at R = 775 when A o=.002, 
and at R = 635 when A o=.0066 (which is well before the onset of the primary wave), reaching 
maximum amplification factors of 8, 15, and 30 respectively. Note that when the initial ampli- 
tude A 0 is large enough, the induced instability can be so strong and secondary instability 
occurs directly by-passing the usual growth of the primary wave. In a situation like this, tran- 
sition prediction schemes based on linear primary theory fail completely. 
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In boundary-layer flows, the onset of the secondary 3D instability is known to be an 
important feature of the early stages leading to transition. For LFC purposes, one might try to 
completely avoid or delay this instability by using suction. Then, one faces questions like 
where to apply suction, and whether it should start before or after the onset of the secondary 
instability, and what is the optimum suction needed to keep the flow laminar. To answer these 
questions, calculations were made to compare the stability characteristics of different cases of 
Blasius flow with continuously applied suction starting at five different locations. These results 
are reported in Fig (5) which shows the variation with R of the growth rates of both the pri- 
mary and secondary disturbances at fixed initial amplitude of the primary wave (Aq=0.001). 
Calculations show that in order to control the secondary 3D instability, suction should be 
applied further upstream near R, of the primary wave and not to apply it near the onset of the 
secondary instability. While investigating the effect of suction on primary TS waves, Reed and 
Nayfeh [27] and Saric and Reed [28] reached similar conclusions that suction should be con- 
centrated not in the region of maximum growth but near its first neutral stability. 

Figure (6) shows the variation with suction parameter r 0 of the maximum growth rates 
with respect to R of both the primary wave and secondary subharmonic disturbance for Ao=.001 
and b=. 15. Point P indicates a suction level that completely stabilizes the primary wave, while 
point S indicates a suction level that allows for a limited primary growth but completely stabil- 
izes the secondary disturbance ( point S is extrapolated due to doubtful Floquet theory results 
when the secondary growth rates are small). An optimum suction requirements for an LFC 
system may be located somewhere upstream of point S allowing for a considerable growth of 
the most unstable primary wave, as well as a limited growth of the secondary disturbance as 
long as the primary amplitude is small and the induced secondary disturbance is not strong 
enough for nonlinear self- and cross-interaction with the primary wave. As A 0 increases, point 
S is expected to move towards point P and cross it over for high enough primary initial ampli- 
tudes. In such situation, suction level required to fully stabilize the primary wave may not 
control the amplification of the secondary disturbance, and very high suction levels are needed 
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to control the flow. 

4.2. Falkner-Skan Flow Controlled by Suction 

The need to control the secondary instability in a boundary layer with pressure gradient is 
more realistic such as the case of flow on airfoils. Herbert and Bertolotti [29] studied the 
secondary instability of Falkner-Skan flows and found that favorable pressure gradient limits 
the band of unstable spanwise wavenumbers, while small adverse pressure gradient is strongly 
destabilizing. This means that larger suction rates are needed to control the boundary layer in 
these cases. For a boundary layer with pressure gradient and suction, the variation of the 
growth rates of the secondary 3D subharmonic with the spanwise wavenumber parameter b 
exhibits the same features given before in Fig(l). 

Figure (7) gives an overall view of the effect of both pressure gradient and suction param- 
eters on the secondary 3D subharmonic at fixed A 0 using the maximum amplification factor as 
a basis for comparison. As adverse pressure gradient increases, the secondary subharmonic 
disturbance becomes more unstable. With the suction parameter increasing, the curves con- 
verge rapidly to lower amplification factors indicating that the secondary instability is more 
sensitive to suction as adverse pressure gradient increases. Figure (7) shows also the increase 
in the maximum amplification factors for different values of A 0 at P<p0. The sensitivity of the 
secondary instability to suction appears to be higher with the increase in adverse pressure gra- 
dient than the increase in A 0 . 

4.3. Mean Profile, Mode Shape, and Vortical Structure 

Modifications to the Blasius profile due to suction and adverse pressure gradient are 
shown in Fig (8). Suction leads to a fuller mean [/-velocity profile and decrease in the magni- 
tude of the V' -velocity. While adverse pressure gradient makes the [/-velocity profile more 
inflectional. At fixed amplitude of the primary wave, these modifications have moderate effect 
on the secondary instability, see Fig (2) for the case of suction, and Ref. 27 for the case of 
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pressure gradient. 

Figure (9) gives a comparison of the mode shapes of the v-, and w -velocity com- 
ponents of the secondary disturbance at different suction and pressure gradient parameters. 
The v -component is small and is slightly affected by suction or pressure gradient. The 
corresponding mode shapes of the u- and v -velocity components of the primary wave are 
shown in Fig (10) for comparison. Suction tends to move the critical layer closer to the wall 
with the maximum of the primary wave velocity components and the maximum of the secon- 
dary disturbance velocity components following it. The location of the critical layer on the t|- 
axis is marked in both figures. By increasing suction, the thickness of the boundary layer 
decreases, and the disturbance is confined to a region closer to the wall where dissipative 
action is strong, thereby increasing the stability of the flow. As adverse pressure gradient 
increases, the opposite occurs and a change over from viscous to inflectional instability takes 
place [29]. 

A more detailed description of the physics of the process of suction control can be 
obtained from contour plots of the vorticity components. Figures (11)-(14) give contours of 
vorticity components in a wave-fixed coordinate system for different cases of suction and pres- 
sure gradients at =1050. In each figure, frames (a)-(d) show the effect of the mean-flow 
modification on the vorticity contours (4=0.04 and B=1 are fixed), while frame (e) shows the 
total effect of the suction on the vorticity contours due to modifications in the mean flow and 
the reduction in the primary amplitude (4=0.006 and B =0.005). These values of 4 and B are cal- 
culated at /? =1050 assuming 4 0 =0.001 at the onset of the primary wave, and B^lxlQr 6 at the 
onset of the secondary subharmonic. Comparison of frames (a) and (d) in each figure shows 
that suction with r<j=-.2 applied to an inflectional profile with Po=.l9 qualitatively produces the 
same vorticity contour plots as for a blasius profile. 

An array of streamwise-periodic concentrations of vorticity is established by the primary 
viscous instability. The strength of the vorticity intensifies with the increase of the amplitude of 
the primary wave. Figure (11) shows plots of the initial 2D vorticity contours of the basic 
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flow ( only spanwise component exist), which peaks near the wall, and extends to the critical 
layer. The location of the critical layer is defined by tick marks on the T|-axis. 

With the onset of the secondary instability, 3D vortical structure is induced by the defor- 
mation of the initial 2D vorticity. Figure (12) shows plots of the spanwise vorticity C, z -contours 
in the x-y plane at z= 0 for different cases of suction and pressure gradient. The figure that is 
plotted over four primary wavelengths (X x =2n/a), indicates that regions of concentrated span- 
wise vorticity are convected downstream, pulled towards the wall as suction increases, and 
lifted away from the wall as adverse pressure gradient increases. In both situations, the con- 
centrated vorticity follows the critical layer as the suction and/or pressure gradient changes. 

The spanwise-velocity variations produce streamwise vorticity C, that is contoured in Fig 
(13). The figure is plotted in the z-y plane at jc=0 over two spanwise wavelengths of the 
secondary disturbance (X*= 27 t/p). It shows an array of counter-rotating vortices extending away 
from the wall as adverse pressure gradient increases, and pulled towards the wall as suction 
increases. 

Following the onset of the secondary instability, spanwise and streamwise vortices experi- 
ence a process of stretching and tilting as they move downstream. Together, they form a 
large scale 3D structure ( A-shaped) which was observed experimentally and numerically. For 
different suction and pressure gradient parameters, Fig (14) shows the deformation in the total 
vorticity (VC*+C?) in an x-z plane almost at the relocation of the critical layer. 

5. CONCLUDING REMARKS 

Calculations show that the secondary 3D subharmonic instability is very sensitive to and 
can be controlled by suction. Meaning, the onset of the instability is delayed, the growth rates 
and amplification factors are reduced, the unstable band of spanwise wavelengths is limited, 
and the vortical structure is closer to the wall. As adverse pressure gradient increases, the sen- 
sitivity to suction increases. For higher initial amplitude of the primary wave, suction is less 
effective in controlling the secondary 3D subharmonic instability. 
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Calculations show that elimination of the primary wave by suction at an arbitrary stage of 
its growth might not delay or eliminate the secondary instability specially if it has already 
taken place. The suction need be applied near R, of the primary wave to gain more control on 
the secondary growth. An optimum suction amount may be lower or higher than the suction 
needed to fully eliminate the primary wave. This is very much dependent on the initial ampli- 
tude of the primary wave. 

Application of suction influences both the mean profile and the primary wave. The effect 
of a modified mean profile on the secondary 3D instability is moderate, while the main effect 
is due to enhancing the growth of the primary wave. By applying suction and including both 
effects, calculations indicate that the most unstable secondary subharmonic disturbance occurs 
at lower spanwise wavenumber. Then, at certain downstream location, an observed flow struc- 
ture of the H-type might be altered to C-type with larger spanwise wavelength or might com- 
pletely disappear as suction increases. 

Evaluation of the effect of suction on the subharmonic secondary instability is a step 
towards the goal of optimizing an LFC system. The idea of monotoring the 3D A-shaped 
structure as its spanwise wavelength changes with suction may be used for that purpose. 
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APPENDIX : Nonzero Elements of the Coefficient Matrices a**, b m , and c m of Eqs (18) 
and (21) : 

ai2=l, a 2 i = ii?(at/-co)+a 2 , a 23 = RDU, a 2A = iaR, 

a 3i = -i a , a 42 = -i a JR , a 43 = -i (at/-io)-a 2 /l? 


b 23 = RDU 


b u =(y+^iar)R 


* 12 = 1 • 


*21 — 1 » 
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z>31 - -b-ifjR , ^34 - - P * 645 = 1 

654 = (V--^iar)U-^i<o)R - (v-^ia r ) 2 + p 2 


b 56 — —R P , b 62, - bn/R , £>63 = -b^R , b $ 5 = — p//? , 


C21 = . 


c 22 = Rvi , 


C 23 - R Du x ,C 54 — (Y ^ I Ctr)/?!* ! , 


Css = £? Vi , 


C 61 = (Y-y*a r )vi , 


C63 = — C 54 /R—DV 1 , C64 = PVl 
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Fig. I Effect of the amplitude A of the primary wave on the growth rate of the 
secondary 3D subharmonic at different suction parameters, R = 1050, and F= 60. 



b 


Fig. 2 Effect of mean-flow modifications on the growth rate of the secondary 3D 
subharmonic at R = 1050 ( curves a, b, and c, with fixed A = 0.01) and the total effect of 
suction (on mean flow and primary wave, curve d , with A - 0.0015). A 0 = 0.00025, and 
F = 60. 
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Fig. 3 Effect of suction parameter on the onset and growth rates of the 
primary wave an( j secondary 3D subharmonic, A 0 = 001, b=.l5, and F-60. 
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In B/Bq 



Fig. 4 Effect of the initial amplitude A 0 of the primary wave on the onset and 
amplification factor of the secondary 3D subharmonic, dotted arrow: onset of the pri- 
mary, solid arrow: onset of the secondary subharmonic, ro=-.05, b=. 15, and F=60. 


-24- 


growth 
rate, 
X 103 



R 


Fig. 5 Continuous suction starting at different locations and its effect on the onset and 
growth rate of both the primary wave and the secondary 3D subharmonic, ^qf.OOI, 
!>=.15, and constant suction level (corresponding to r 0 =-.05 at /t = l000). 


-25- 


Max 
growth 
rate, 
X 103 



Fig. 6 Variation of maximum growth rate with the suction parameter for both the pri- 
mary wave and the secondary 3D subharmonic, /to= 001, b = . 15, and F=60. 



Fig. 7 Effect of suction parameter on the maximum amplification factors of both the 
primary wave and the secondary 3D subharmonic at various pressure gradient parame- 
ters, /tp=.001. b=A5, and 5=60. 
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Fig. 8 Effect of suction and pressure gradient parameters on mean flow profiles. 


- 27 - 



Fig. 9 Effect of suction and pressure gradient parameters on the eigenmodes of the 
secondary subharmonic disturbance, * = 1050, A = .0l, b=. 15, and F= 60. 



Fig. 10 Effect of suction and pressure gradient parameters on the eigenmodes of the 
primary wave, R = 1050 and F= 60 . 
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(a) r 0 =O, po=0, A = .04 


(b) r 0 =-.6, 3o=0, A = .04 
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Fig. 11 Initial 2D vorticity contours of the basic flow at different suction and pressure 
gradient parameters. Flow direction is from left to right, R = 1050, F=60. 
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(C) r 0 =0, p 0 =-.19, A = .04, 5=1 


(d) r 0 =-.2, Pq =-.19, A = .04, 5 = 1 


(e) r 0 =-.05, Po =0, A = .006, 5 = .005 



Fig. 12 Spanwise vorticity contours of the 3D flow in the x-y plane at z=0 for different 
suction and pressure gradient parameters. R = 1050, b=. 15, and F=60. 


- 31 - 



(a) r 0 =0, P 0 =0, A = .04, 3 = 1 


(b) r 0 =-.6, Po =0, A = . 04, 3 = 1 



(C) r 0 =0, Po =-.19, A = .04, B=l 
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(e) r o =-.05, Po =0, A = .006, 3 = .005 
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Fig. 13 Streamwise vorticity contours of the 3D flow in the z-y plane at x=0 for 
different suction and pressure gradient parameters. 3 = 1050, b=. 15, and /= 60. 
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(a) r 0 =0, po =0, A = .04, B = 1 (b) r 0 =-.6, Po=0, A = .04, fl = 1 



Fig. 14 Total vorticity contours of the 3D flow in the x-z plane almost at the y- 
location of the critical layer for different suction and pressure gradient parameters. 
R = 1050, b=. 15, and ^=60. 
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